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Abstract 

■ A new definition of Lie invariance for nonlinear multi-dimensional boundary value 
. problems (BVPs) is proposed by the generalization of known definitions to much wider 
I classes of BVPs. The class of (l-|-3)-dimensional nonlinear BVPs of the Stefan type, 

^ ■ modeling the process of melting and evaporation of metals, is studied in detail. Using 

O . the definition proposed, the group classification problem for this class of BVPs is solved 

^ I and some reductions (with physical meaning) to BVPs of lower dimensionality are made. 

■ Examples of how to construct exact solutions of the (l-|-3)-dimensional nonlinear BVP 
^ ■ with the correctly-specified coefficients are presented. 



1 Introduction 



Currently, Lie symmetries are widely applied to study partial differential equations (PDEs) 
(including multi-component systems of mult i- dimensional PDEs), notably, for their reductions 
to ordinary differential equations (ODEs) and constructing exact solutions. There are a vast 
>- 1 number of papers and many excellent books (see, e.g., [T145] and references cited therein) devoted 

^ ■ to such applications. However, one may note that the authors usually do not pay any attention 

1/-^ ■ to the application of Lie symmetries for solving boundary value problems (BVPs). To the best 

! of our knowledge, the first papers that did so were published at the beginning of 1970s (see [B] 

and |7] and their extended versions presented in books [S] and [2], respectively). The books, 
^ I which highlight essential role of Lie symmetries in solving BVPs and present several examples, 

were published much later [2ll9|[T0]. 

The main object of this paper is a class of (1-1-3) -dimensional nonlinear BVPs of the Stefan 
^ '. type. These problems are widely used in mathematical modeling of a wide range of processes, 

I which arise in physics, biology, chemistry and industry [TTHTB] . Nevertheless, these processes 

can be very different from the formal point of view, they have the common peculiarity, unknown 
moving boundaries (free boundaries). Movement of unknown boundaries is described by the 
famous Stefan boundary conditions [T614T8] . It is well known that exact solutions of BVPs 
of the Stefan type can be derived only in exceptional cases and the relevant list is not very 
long at the present time (see [TT |[TUH25] and references cited therein). Notably, those exact 
solutions were constructed under additional conditions on their form and/or the coefficients 
arising in the relevant BVP. It should also be stressed that all analytical results derived in those 
papers are concerned with two-dimensional BVPs. To the best of our knowledge, there are no 
invariant solutions (with physical meaning) of multidimensional BVPs with free boundaries 
excepting particular problems with radial symmetry, for which analytical solutions are found 
(see, e.g., [TTIET] and references cited therein). Perhaps |28j (see also references by the same 
author cited therein) is unique because it contains such exact solutions for (l-|-3)-dimensional 
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hydro dynamical problems with free boundaries. Of course, there are many interesting papers, 
devoted to the rigorous asymptotic analysis of such BVPs, leading to the relevant analytical 
results (see, e.g., [29ti3T] and references cited therein). 

From the mathematical point of view, BVPs with free boundaries are more complicated 
objects than the standard BVPs with fixed boundaries. In a particular case, each BVP with 
Stefan boundary conditions is nonlinear; nevertheless, the basic equations may be linear p^l32] . 
Thus, the classical methods of solving linear BVPs (the Fourier method, the Laplace transfor- 
mation and so forth) cannot be directly applied for solving any BVP with free boundaries. On 
the other hand, it can be noted that the Lie symmetry method can be more applicable just 
for solving problems with moving boundaries than for other BVPs. In fact, the structure of 
unknown boundaries may depend on invariant variable(s) and this provides the possibility to 
reduce the BVP in question to that of lover dimensionality. This is the reason why some authors 
have applied the Lie symmetry method to BVPs with free boundaries [6|[7tl2 H [28 |l33] - [35] . 

The paper is organized as follows. In section 2, we propose a new definition of Lie invariance 
for any BVPs with basic evolution equations, which generalize the known definitions, and 
formulate the algorithm for solving the group classification problem. In section 3, we apply 
the definition and the algorithm to the class of (l+3)-dimensional BVPs of the Stefan type, 
used to describe melting and evaporation of materials in the case when their surface is exposed 
to a powerful flux of energy. The main result is presented in Theorem 2, which is a highly 
non-trivial generalization of that, derived for two-dimensional BVPs in [23] . In section 4, 
all possible systems of subalgebras (optimal systems of subalgebras) for a subclass of (1+3)- 
dimensional BVPs admitting a five-dimensional Lie algebra of invariance are constructed. We 
reduce such problems to the two-dimensional BVPs via the non-conjugate two-dimensional 
subalgebra. Moreover, we show that this reduction admits a clear physical interpretation. 
Examples of how to construct exact solutions of the (l+3)-dimensional nonlinear BVP with 
the specified coefficients of diffusion are presented too. Finally, we present conclusions in 
section 5. 



2 Definition of Lie invariance for a BVP with free bound- 
aries 

Consider a BVP for a system of n evolution equations [n > 2) with m+ 1 independent {t, x) 
(hereafter x = {xi,X2, ■ ■ ■ , Xm)) and n dependent u = {ui,U2, ■ ■ ■ , Un) variables. Let us assume 
that the basic equations possess the form 

= (x, u, 4^')) , z = 1, . . . , n (1) 

and are defined on a domain (0, +oo) x Q G M™+^, where Q is an open domain with smooth 
boundaries. Hereafter, the lower subscripts t and x denote differentiation with respect to these 
variables and u^''' denotes a totality of partial derivatives with respect to x of order ki (for 
example, ^ = m^,, . . . , Ua;^). 

Consider three types of boundary and initial conditions, which widely arise in applications: 

Sait,x)=0: Bi(t,x,u,Ua:,...,ui^''''^^ =0, a = l,...,p, j = l,...,na, (2) 
Shit,x) = 0: Bi{t,x,u,...,uf''\si'\...,sf''^)=0, b=l,...,q,l = l,...,nh, (3) 
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and 

7c(t, x) = oo : [t, X, u,Ux,. ■ u^x""^) = 0, c = 1, . . . , r, m = 1, . . . , ric- (4) 

Here < k = max{A;i, . . . , A;„}, k[ < k, k^ < k and Kl are the given numbers, Sa{t, x) and 
7c(t,x) are the known functions, while the functions Sb{t,x) defining free boundary surfaces 

must be found. In ([3]), the function S^^ '' denotes a totahty of derivatives with respect to t and 
X of order (for instance, S^^^ = ||^, . . . , §§^)- We also assume that all functions arising 
in are sufficiently smooth so that a classical solution exists for this BVP. 

We note that boundary conditions (jlj) essentially differ from those ([2]) because they are 
defined on the non-regular manifolds Aic = {icitjx) = oo}. Such conditions appears if one 
considers EVPs in the unbounded domains and often leads to difficulties. For example, one may 
check that the definition of BVP invariance presented in [Tl|2] is not valid for such conditions. 

Consider an A^-parameter (local) Lie group Ga^ of point transformations of variables {t, x, u) 
in the Euclidean space (open subset of defined by the equations 

t* = T{t,x,e), X* = Xi{t,x,e), u* = Uj{t,x,u,e), i = 1, . . . ,m, j = 1, . . . ,n, (5) 

where e = {ei,e2, ■ ■ ■ ,eiy) are the group parameters. According to the general Lie group theory, 
one may construct the corresponding A^- dimensional Lie algebra L^r with the basic generators 



dT{t,x,e) 



1 So 



dXi{t,x,e) 



j _ dUjit,x,u,e) 

e=0 ' 



where f " = ^ ' . - 

In the extended space of the variables {t, x, u, S) (hereafter S = (Si, Sg) ), the 

Lie algebra defines the Lie group Gn- 

t* = T{t,x,e), X* = Xi{t,x,e), u* = Uj{t,x,u,e), = Sb, (7) 

where i = 1, . . . ,m, j = 1, . . . ,n, b = 1, . . . ,q. 

Now we propose a new definition, which is based on the standard definition of differential 
equation invariance as an invariant manifold Ai and generalizes the previous definitions of 
BVP invariance (see, e.g., [Il|2l[33l|36]). 

Definition 1 The BVP of the form is called invariant with respect to the Lie group 

Gn 0^f 

(a) the manifold determined by equation (QP in the space of variables (^t, x,u, . . . ,ul^^^ is 
invariant with respect to the kth- order prolongation of the group Gn', 

(b) each manifold determined by conditions ^ with any fixed number a is invariant 
with respect to the kath-order prolongation of the group Gn in the space of variables 

{t, x,u,. . ., ui^"^^ , where ka = max{kl, j = 1, . . . , Ua}; 

(c) each manifold determined by conditions ^ with any fixed number b is invariant 
with respect to the kbth-order prolongation of the group Gn in the space of variables 
(t, x,u,..., ui^''\ Sb,..., Sj^''''^ ) , where kb = max{/c[, Kl, 1 = 1,..., nj; 
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(d) each manifold determined by conditions ^ with any fixed number c is invariant 
with respect to the kcth-order prolongation of the group Gn in the space of variables 

(^t,x,u, . . . jU^x"^^ , w/iere A;c = max{/c™, m = 1, . . . , ric}. 

Definition 2 The functions uj = x), j = 1, . . . ,n and Sb = x), b = 1, . . . ,q form an 
invariant solution (w, S) of the BVP of the form (CP^® corresponding to the Lie group ^ if 

(i) {u,S) satisfies equations (dj and conditions (CP-Q'; 

(a) the manifold Ai = {uj = ^j{t,x), j = Sb = b = l,...,q} is an 

invariant manifold of the Lie group Ij^. 

Remark 1 Definition 1 can be straightforwardly generalized on EVPs with governing systems 
of equations of hyperbolic, elliptic and mixed types. However, one should additionally assume 
that n-component governing system of PDEs is presented in a 'canonical' form (some authors 
uses the natation 'involution form' in this context), i.e. one possesses the simplest form and 
there are no non-trivial differential consequences. 



If the system of differential equations contain arbitrary functions as coefficients (formally 
speaking, they can be constants), the group classification problem arises. Such problems was 
formulated and solved for a class of non-linear heat equations (NHEs) in a pioneering work [37] 
(see also [5]). At the present time, there are algorithms for rigorous solving group classification 
problems (see, e.g., |3S] and references cited therein), which were successfully applied to different 
classes of PDEs. Thus, if system ([T]) and/or the boundary conditions contain arbitrary 

functions as coefficients, then we should formulate and solve the group classification problem 
for the BVP of the form 

We propose the following algorithm of the group classification for the class of BVPs 

(I) to construct the equivalence group i?eq of local transformations, which transform the 
governing system of equations into itself; 

(II) to extend the space of i^eq action on the variables S = {Si, Sq) by adding the identity 
transformations for them (see the last formula in (^^) and to denote the group obtained 

as L^G<\) 

(III) to find the equivalence group E^^^ of local transformations, which transform the class of 

BVPs into itself, one extends the space of the Eeq action on the prolonged space, 

where all arbitrary elements arising in boundary conditions ([2])-(jl]) are treated as new 
variables. 

(IV) to perform the group classification of the governing system ([1]) up to local transformations 
generated by the group Ef^^; 

(V) using Definition 1, to find the principal group of invar iance G^, which is admitted by each 
BVP belonging to the class in question; 

(VI) using Definition 1 and the results obtained at step (IV), to describe all possible Ef^^- 
inequivalent BVPs of the form admitting maximal invariance groups of higher 

dimensionality than G'^. 
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Figure 1: A scheme for the process of melting and evaporation of a sohd material which is 
exposed to a powerful energy flux. 

3 Lie invariance of a class of (l+3)-dimensional nonlin- 
ear BVPs with free boundaries 

3.1 Mathematical model of melting and evaporation under a pow- 
erful flux of energy 

Consider the process of melting and evaporation in half-space = {x = (xi, X2, X3) : X3 > 0} 
occupied by a solid material, when its surface (initially it is the plane X3 = 0) is exposed to a 
powerful flux of energy. We neglect the initial short-time non-equilibrium stage of the process 
and consider the process at the stage when three phases already take place and assume that 
this occurs at any moment t G T = (t^,, -|-cxo), where t^, is a positive real number. Thus, the 
heat transfer domain Q{t) = f2 x T consists of three sub-domains occupied by the gas, liquid 
and solid phases, which will be denoted by ^^o(^), and ^2(^)1 respectively, and the phase 

division boundary surfaces, S'i(t,x) = and S'2(t,x) = (see Fig.l). In other words, the 
domain Q{t) admits the disjoint decomposition 

Q{t) = Qo{t) u ri(t) u Qi{t) u r2(t) u Q2{t), 

where 

Tkit) = {(t,x) : 5fc(t,x) = 0, t G T, X G fi}, k = l,2, 

no{t) = {(t,x) : 5i(t,x) < 0, S2{t,x) < 0, t G T, X G fi}, 

ni{t) = {(t,x) : 5i(t,x) > 0, S2it,x.) < 0, tel, X G fi}, 

^2(1) = {(t,x) : 5i(t,x) > 0, S2it,x.) > 0, t G T, X G r^}. 

Let us consider a class of (l-|-3)-dimensional nonlinear BVPs of the Stefan type used to 
describe melting and evaporation of materials in the case when their surface is exposed to a 
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powerful flux of energy p ^ [T5 | fT8 | l39]: 



V(Ai(ri)VTi) = Ci(Ti)^, (t,x)Gr]i(t), (8) 
V{X2{T,)VT2) = C,iT,)-^, (t,x)Gfi2(t), (9) 

5i(t,x) = : Ai(r,)-^ = H,v, ■ m - Q(t) ■ m, t, = t,, (lo) 

aril 

BT BT 

S2it,X.) = : X2iT„,)-^ = Ai(r„)-^ + HrnV2 " Ha, Ti = T2 = T^, (11) 

|x| = +00 : T2 = Too, t G 1, (12) 

where Ty, Tm and Too are the known temperatures of evaporation, melting and solid phases of the 
material, respectively; Xk(Tk), k = 1,2 are the positive thermal conductivities; Ck(Tk), H^, Hm 
are the positive specific heat values per unit volume; Q(t) = {Qi{t),Q2{t),Q3(t)) is the energy 
flux being absorbed by the material; Sk(t,x.) = 0, A; = 1,2 are the phase division boundary 
surfaces to be found; Vfc(t, x), A; = 1, 2 are the phase division boundary velocities; n^, A; = 1, 2 
are the unit outward normals to the surfaces Sk{t,x.) = 0, A; = 1,2; Tfc(t,x), k = 1,2 are the 

unknown temperature fields; V = ^^f^, ^fj, gf^j i the subscripts k = 1 and k = 2 correspond 
to the liquid and sohd phases, respectively. 

Here equations ([8]) and ([9]) are basic and describe the heat transfer process in liquid and 
solid phases, the boundary conditions (fTOj) present evaporation dynamics on the surface 5*1 = 0, 
and the boundary conditions (fTTl) are the well-known Stefan conditions on the surface 6*2 = 
dividing the liquid and solid phases. Since the liquid phase thickness is considerably less than 
the sohd phase thickness, one may use the Dirichlet condition ( TT^ . where Too can be treated 
as the initial temperature of material. We also assume that the gas phase does not interact 
with liquid and solid phases; hence the problem in question does not involve any equation for 
the gas phase. 

From the mathematical and physical points of view, we should impose the same additional 
conditions on the functions and constants arising in the EVP class in question, which guarantee 
existing classical solutions. Namely, we assume that all functions in (l8|)- (fT2|) are sufficiently 
smooth; the free boundary surfaces Sk{t,x) = satisfy the restrictions ^ 7^ and |VS'fc| 7^ 
0, k = 1,2; the projection of the heat flux vector Q(t) on the normal ni is nonzero, i.e. 
Q(t) • 111 7^ 0, and Vfc ■ 7^ 0, A; = 1,2. Finally, the constants Ty, Tm and Too satisfy the 
natural inequahties Ty > T^ > Too. 

First of all, we simplify the governing system of equations ([8]) and ([9]) using the Goodman 
substitution 

Ti T2 

u = 01 (TO = / Ci(C) dC, V = MT2) ^ ! ^2(0 di. (13) 



Substituting (|T3|) into (l8|)-( fT2|) and making the relevant calculations, we arrive at the equivalent 
class of EVPs 

Bii 

_ = v(di(«)Vw), (^,x)G^]l(t), (14) 

^^=V{^d2{v)Vv), {t,^)en2{t), (15) 
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du 

Si{t, x) = : diy- — = H^Vi ■ ni - Q(t) ■ rii, u = Uy, (16) 

OH I 

dv du 

S2{t, X) = : d2m-g^ = ^Im^^ + -^^^2 ■ n2, U = Um, V = Vm, (17) 

|x| = +00 : V = Voo, t G T, (18) 

TjYi Tm Too 

where Uy = J Ci(C)rfC, Mm = / Ci(C)rfC, Vm = J C2{0d^, = J 6*2(0 M^) = 



Ci(^'-'(I)) ' '^^(v) = ^^^4^' '^^'^ ^ di{uy), dim = di{um), d2m = d2{vm) (here A; = 1,2 are 
the inverse functions of 0^, the functions di{u) and d2{v) are strictly positive and ^ Um, 



-'OO J 



3.2 Group classification of the basic equations (1141)— (1151) 

One sees that the BVP of the form f[T^ - f[TH]) consists of two standard NHEs with arbitrary 
smooth functions di{u) and d2{v) and the boundary conditions ffTB]) - ffT5]) containing an arbitrary 
vector function Q(t) and a number of arbitrary parameters. Thus, we deal with a class of BVPs, 
and to carry out the group classification the algorithm formulated in Section 2 can be used. 

According to item (I), we need to find the group of equivalent transformations of the non- 
coupled system of NHEs f lT^ -( lTB]) . Note that this group is well-known in the case of a single 
NHE (see, e.g., [12] )• However, one cannot extend this result in the case of system f lT^ - fllSp 
in a formal way because the group obtained may be incomplete. Thus, we carefully check the 
group of equivalent transformations for the non-coupled system f lT4|) -f fT5|) . 

Lemma 1 The equivalence transformations group Ecq of system [I4\)- I(T5\} consists of the group 
of continuous equivalence transformations 

t = at + 7o, X = (3A,ii3i)A^ii32)AM)X + T (i, J, k = 1,2,3; i =^ j,i ^ k, j ^ k), 

- 

u = 6iu + 74, V = 62V + 75, di = — di, d2 = — d2, 

a a 

and the group of discrete equivalence transformations 

1) t ^ —t, Xj — 7- (— l)-'Xj (z = 1, 2, 3, j = 0, 1), u —i- —u, t; — )■ — 1>, di^—di, ^2 ~> — <^2;' 

2) t ^ t, — 7- {—lyxi {i = 1, 2, 3, j = 0, 1), u ^ V, t> — )■ u, di — )■ d2, d2 ^ di. 

Here a > 0, /3 > 0, . . . , /^a, 70, . . . , 75, 5i > 0, ^2 > are arbitrary constants; 

/ cos^ sin^ \ / COS0 sin0 \ / 1 \ 

Ai{e)=\ -sine cose ,A2(^)= 1 ,^3(0)= cos^ sine 

\ 1 / \ - sin ^ COS / \ - sin 6* cos ^ / 

are the matrix of rotation in space, and 





f ^1 ) 




r7i\ 




X2 




72 








^73/ 
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Proof. To find the group £^eq, we use the standard infinitesimal method [3113], i.e., search 
for the generators of the form 



Y 



d 



d 



d 



+ TJ^ h h yU"^ 

dxa du dv ddi 



where a = 1,2,3 (hereafter, summation is assumed from 1 to 3 over repeated indices a or b). 
The generator Y defines the group Seq of equivalence transformations 

t = 0(t,X, = Va(^,X, li,!)), M = $1 (t , X, M, t;) , V = <l>2 (t, X, M, i)) , 

di = \E'i(t,x, -u,-u,di,(i2), d2 = '^2{t,^,u,v,di,d2) 
for the class of systems (|T^ - (|T5|) iff Y obeys the condition of invariance of the following system: 



^-V{d^Vu) 
|-V(4V.) 



0, 
0, 



ddi 
'dt 

dd2 



0, 



0, 



ddi 

dXa 

dd2 



0, 



0, 



ddi 
dv 
dd2 



0, 



0. 



(19) 



(20) 



dt dxa du 

Here, the variables u and v are considered in the space (t,x), while dk in the extended space 
(t, X, M,f). The coordinates and rf of the operator Y are the functions of the variables 

t,x, w,!), while the coordinates /i'^ are the functions of t,:K,u,v,di,d2- Thus, the invariance 
criterium for system fll9p -( 12II]) is given by the formulae 

'ddi 



0, 



y(2) /'^ 
dt 



V {d2Vv) 



0, r(^) 



dt 

dd2 
~dt 



= 0, ( ) =0, y^^^ I ) =0, 

(21) 

\UXa/ g \'-'^ / S 

(22) 

where is the manifold, defined by system ( IT9|) - (l20|) . y*^^) is the second prolongation of the 
operator Y calculated via the known formulae |5lll3]. Using these formulae and ( l2T]) - (l22|) and 
carrying out the relevant calculations one obtains the 13-dimensional Lie algebra with the basic 
operators 

Yi = ty. = ^..y^ ' ' 



Y. 



X2 



d 

dxi 



- Xi 



Y 
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df 
d 
dx2 
d_ 
du' 

t—-d — 
dt ^ ddi 



dxi 



— d'}- 



dv ' 
d 



d 

dxi 

Yi, 



d 



- Xi 



u 



dxs 
d_ 
du' 
d 



, Yr 
Yii 



X3 



d 
dx2 
d_ 
dv' 



- X2 



_d_ 

<9xa 



+ 2di— + 2d2 



ddi 



ddo 



dd2 dXa 

One easily checks that this algebra generates the group Seq- 

Finally, to obtain the group of equivalence transformations E^q, we should add the discrete 
transformations 1) and 2) listed in Lemma 1. It is easy to verify by the direct calculations that 
system (|T4|) - (|T5|) is invariant under these discrete transformations. 

The proof is now complete. ■ 

According to item (II) of the algorithm, we should now construct the group E^q by adding 
the identity transformations 5*1 = <S'i,5'2 = <S'2. To realize item (III), one needs to apply the 
transformations generated by the group Eeq to the boundary conditions (|T6|) - (|T8l) and to find 
the group Ef^^. The result obtained is formulated as follows. 
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Table 1: Lie algebras of the NHE system ( [I^ - (|T5|) {ki,k2,m and n are arbitrary non-zero 
constants, 6 < a = 1,2,3; while a{t,x) and j3{t,x) are arbitrary solutions of the linear heat 
equations at = kia^x and Pt = k2(5xxi respectively.) 



no 


di{u) 


d2{v) 


Basic operators of MAI 




1. 


V 


V 


AE{\,'i) = {dt, dx,,Xadx^ - xtdx^,2tdt + Xad^,) 




2. 


ki 


V 


AE{l,3),udu,a{t,x)du 




3. 


V 




AE{l,3),vd,,/3it,^)d, 




4. 


e" 




AE{l,3),Xadx^ + 2du + 2d, 




5. 






AE(l,3),Xnd^ +2d„ + -vd„ 




6. 




e" 


AE{l,3),Xadx^ + ^udu + 2d, 




7. 






AE{1, 3),D = Xadx^ + ^udu + ivd. 




8. 


4 

M 5 


4 

V 5 


AE{1, 3), xpc?^^ — 2xbXadx^ + 5xbudu + dx^vd^, D with m = n = - 


4 
"5 


9. 






AE{l,3),udu,vdy,a{t,:s.)du, f3(t,:s.)d^,,Ga = td^, - Xa {^ud^ + 

n = t^dt + tXadx^ - ^i^dxp + ^k^t)udu - 4|;(|X|2 + ^k2t)vd. 




10. 


ki 




AE{\, 3),udu, vdv, vdu, ud^, a{t, x)5„, (3(t, x)9^, Ga and 11 with ^2 = 


h 



Lemma 2 The class of B VPs admits the group of equivalence transformations Ef^^ : 

i=at + jo, X = (3A{(3,)Aj{P2)Ak{(3s)X + T {i,j,k = 1,2,3; i ^ j,i ^ k, j ^ k), 

~ ^2 ^ ^2 

u = Siu + 74, V = S2V + 75, 5"! = 5*1, 5*2 = 5*2, di = — di, d2 = — d2, 

a a 

rl -1 ri ^ -Aw n --n n --H 

'^Iti r- ^Iti) "^Im e ""Irrt) W2m r- '^2mi n vi n ^i 

Oi Oi O2 p p 

= diU^ + 74, Ura = ^lUra + 74, Vm = ^2^"™ + 75, ^^oo = S2V00 + 75 

Q = A(/3i)A,(/32)Afc(/33)g iij,k = 1,2,3; i ^ j,i ^ k, 3 ^ k), 
with arbitrary coefficients a, /?, . . . , /^s, 70, . . . , 75, 5i, 82 obeying only the condition 

a^5i52 ^ 0. 

According to item (IV) of the algorithm, we should now perform the group classification of 
the governing system f[T^ - f[T51) up to local transformations generated by the group Ef^^ . The 
result can be presented as follows. 

Theorem 1 All possible maximal algebras of invariance (MAIs) (up to equivalent represen- 
tations generated by transformations from the group Ef^^ ) of system [l4\)- fTS\) for any fixed 
vector function {di, ^2) with strictly positive functions di{u) and d2{v) are presented in TableUl 
Any other system of the form [14^- fT5\) is reduced to one of those with dijjusivities from TableU\ 
by an equivalence transformation from the group E^^^ . 



Proof. First of all, we remind the reader that a complete description of Lie symmetries for 
the class of mult i- dimensional nonlinear reaction-diffusion systems 

— = V idi{u)Vu) + Fi{u,v), 
dv 

— = y[d2{v)Vv) + F2{u,v), 

where Fi and F2 are the arbitrary nonzero smooth functions, was obtained in [mHS] (the 
case of constant diffusivities) and [IH] (the case of non- const ant diffusivities). Nevertheless, the 
detailed examination of the special case Fi = F2 = was omitted in the papers [Hl - HG] . we can 
use the relevant system of the determining equations with Fi = F2 = to find all possible MAIs 
of system (fT4|) - (fT5|) . Let us assume that the MAI in question is generated by the infinitesimal 
operator 

d d d d 

X = f{t, X, u,v)— + C{t, X, u, v)-^ + 77^(t, X, u,v)— + rfit, x, u, v) — , 

where r]^ and 77^ are the unknown smooth functions. The form of the system of deter- 

mining equations to find these functions essentially depends on the derivatives ^ and Let 
us consider the most general case ^ ■ 7^ 0. The relevant system of determining equations 
has the form (see formulae (11)-(16) with Fi = F2 = in [16]) 

=C = C = 0, C = C = vl = vl = 0,a = 1,2,3, (23) 
= 4 = 4' + el = 0, 6 < a = 1, 2, 3, (24) 
e° = 2C -V'^ log d^, = 2C - v'-^ log d2, a = 1,2, 3, (25) 

^P^+^viUogd, = Ae-cndir', 2^+2,i^iogd2 = Ac-cnd.r\ m 



2C-e?-.i-.L(|;log..)" + .^|^log^,a^l,2,3 



(27) 



2C - - = Vl log d,^ " + ^^1; log ^, a = 1, 2, 3, (28) 

r]l - = 0, - Ari^d2 = 0. (29) 

If di{u) and d2{v) are the arbitrary smooth functions, then system (I23|) - (l29|) can be easily 
solved resulting the eight-dimensional Lie algebra y4F(l,3) with the basic operators 

d d d d 

Pt= tt:, Pa= 7^ , Jab = XaPb " XbPa, Dq = 2t— + X^-— , 6 < C = 1, 2, 3. 

at OXa Ot OXa 

Now, we should find all possible pairs of the function di{u) and d2{v) leading to extensions 
of the algebra AE{1, 3). It is evident that equations ( 123|) and (p5l) can be easily solved and the 
functions 

= A{t), C = B''{t,x), a =1,2,3, (30) 
= {2B: - A) (l^logrfi) , = (2i3: - A) (^108^2) , a = 1,2,3 (31) 
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are obtained (here A and 5" are the arbitrary functions at the moment). However, the functions 
5" can be specified using flM|) as follows: 

B"" = coa + CabXb + ka(t)\:>c\'^ -2kb{t)xbXa, a,6 = 1,2,3, (32) 

where Cab + Cba = {a b),cu = . . . = C33 and Cqq, Cab, ka{t), kb{t) are arbitrary functions. 

Substituting formulae fISUl) and fl?Il) in flTTI) and (1251) we arrive at the system of classification 
equations to find pairs of the functions {di, ^2): 

otherwise, the conditions 

2B:-A = 0, a = 1,2, 3, (34) 

must hold. However, substituting ( l32|) in conditions and taking into account (l26l) and 
( 15U]) -( 1?T]1 we immediately obtain only the Lie algebra AE{1,3). System which consists of 
two independent ODEs, can be easy solved: 

'^^^^) = \ Aexp(«in), (35) 



D2iv + C2)"' 
D2exp{aiv), 

where 7^ 0, ^ and are arbitrary constants {k = 1,2). 

Substitutions the functions and r^^ from ( 130|) -( 13T|) in ( 126|) we arrive at the equations 

Since the functions Sq, (a = 1,2,3) depend only on t and x, there are only two possibilities. 
The first one is AS" = and then, applying f l32p . we obtain 

= CabXb + CaO, Cab, CaO E R. 

The second possibility, having AS" 7^ 0, requires 

where mi and m2 are some constants. Using f l35l) - fl36l) it is easily seen that conditions fl37|) can 
be satisfied only for = m2 = 0, and then 

di{u) = Di{u + Ci)--^ , di{v) = D2{v + C2)~~' . (38) 

Thus, we have obtained all possible forms of the functions di{u) and d2{v), namely formulae 
( 15^ -f l5B]) and (EH]) J leading to extensions of the invariance algebra ^^(1,3) of the nonlinear 
system ( 1T4|) -( [T5|) . when both diffusivities are non-constant. 

Finally, taking into account the point transformations from the group Ef^^ , we immediately 
obtain cases 4-8 from Table [TJ Other possibility, when at least one of the functions di{u) and 
d2{v) is constant, was examined in a similar way, and cases 2,3 and 9, 10 from Tabled] were 
obtained. 

The proof is now complete. ■ 

Remark 2 Cases 2 and 5 from TableU\are equivalent to 3 and 6, respectively, if one takes into 
account the discrete transformations 2) from the group E^q. However, these transformations do 
not belong to the Ef^^ , because the boundary conditions ( flgj) and [T^) are not invariant under 
them. 
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Table 2: Lie invariance of BVP (Il4D-([l8]) 



no gi(t) Q2{t) Qsit) MGI 



1. 


V 


V 


V 






2. 








q{t) 


?i,T2,f3,f5 




3. 






Q3 


fi,f2,f3,n 




4. 


(lAiogt) 




13 

Vt 


Ti, fa, fa, TV 




5. 








Q 


To, Ti, T2, T3, 


^5 


6. 








1 


Ti, T2, T3, T4, 


T5 



3.3 Group classification of the class of BVPs (fT4l^-(TT8l) 



We note that each MAI from Table [T] generates the corresponding maximal groups of invari- 
ance (MGI) with the transformations, which can be easily derived from the basic generators 
listed therein. These transformations are well-known and used below. 

Taking this into account and according to items (V) and (VI) of the algorithm, we formulate 
the main theorem giving the complete list of Lie symmetries of the BVP class (fT^ -ljli 



Theorem 2 The BVP of the form [T^-fT8\). with the arbitrary given functions di{u), d2{v) 



(di{u) 7^ d2{v) ) and Qaif), a = 1, 2, 3, is invariant under the three-parameter Lie group (trivial 
Lie group) presented in case 1 of Tahle\^ The MGI of any BVP of the form [T^- ^W\) does not 
depend on the form of di{u) and d2{v). There are only five BVPs from class [l4\)- fT^) with the 
correctly specified functions Qa{t), a = 1,2,3 admitting the MGI of a higher dimensionality, 
namely: four- or five-parameter Lie groups of invariance (up to equivalent representations 
generated by equivalence transformations from the group Ef^^ ). These MGIs and the relevant 
functions Qa{t), a = 1,2,3 are presented in cases 2-6 of Table\M 

Remark 3 In Table\^ the following designations are used: q 0, qs, A are arbitrary constants, 
q{t) ^ is an arbitrary function while the functions 

= Qi COST + q2 sinr, 62 (r) = -qi sinr + q2 cost, 

where qi, q2 are arbitrary constants satisfying the condition ?i + ?| 7^ if r ^ 0. 

The explicit form of the transformations generating the MGI are: 



To 


■.t*=t + eo, 


* * 

X — 1 5 — 5 


xl = X3, 


* 

u 


= u, 


* 

V 


= V, 


SI 




s; 


= S2, 


fi 


• t — X — 


Xi -\- Ei^ X2 — X2-, 


* 


* 

u 


= u, 


* 

V 


= V, 


SI 


= Si, 


S*2 


= S2, 


f2 


. t — ij ^ X-^ — 


^2 = ^2 + ^25 


^3 ~ ^3) 


* 

u 


= u, 


* 

V 


= V, 


SI 


= Su 


S*2 


= 82, 


fs 


. t — i ^ X-^ — 


* * 
X\^ X'2 — ^2 5 — 


X3 + £3, 


* 

u 


= u, 


V* 


= V, 


SI 


= Si, 


s; 


= S2, 


t* 


— ™* _ 


6 1 5 *^2 — ^ '^'2^ 


X* = e'' 


X3, 


u* = 


u, 


V* = 


-- V, 


SI = 


Si, 


S*2 = S, 
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T^-.t* = t, Xl = 61(65), X*2 = 92{£5), X*3 = X3, U* = U, V* = V, SI = Si, S2 = S2, 

Tq : t* = t + £6, = 6i{XeQ), x*2 = 92{XeQ), x*^ = X3, u* = u, v* = v, 5"* = Si, S2 = S2, 
fj-.t* = e^'H, xl = e^'OiiXej), xl = e'^^slAe?), x^ = e''x^, u* = u, v* = v, SI = Si, S; = S2, 
where 

9i (r) = Xl cos T + X2 sin r, (t) = —xi sin t + X2 cos r. 

Proof. Let us consider the case of the arbitrary function di{u) and d2{v). According to 
Theorem 1 (see case 1 of Table 1), in this case the NHE system f[T^ - (IT^ admits the eight- 
parameter group Gg of invariance generated by the groups Tq,Ti,T2,T^,T/i and the rotation 
groups 

Tab ■ t* =t, Xl = -XbSiuEab + Xa COS Eab, xl = Xb COS Sab + XaSiuSab, X*^ = X^, U* = U, V* = V. 

where 6 < a = 1, 2, 3, a, c^h, and Sab are the group parameters. 

Since the BVP of the form f[T^ - f[T51) has two free boundaries, Si{t, x) = and S2{t, x) = 0, 
we need to extend the group G\ by adding identical transformations for the new variables 
S\ = Si and S'2 = S'2. We will denote the obtained group by Gg (the relevant sub-groups will 
be denoted in the same way, for instance, T31). 

By straightforward calculations, one easily checks that the boundary conditions (|T71) and 
f|T8|) are invariant under the group Gg. The situation is essentially different if one examines 
the invariance of the boundary condition (fT6|) with respect to Gg. To simplify calculations, the 
boundary conditions fll6p should be rewritten in the form (see the monograph [13j, P. 18) 

Siit,x) = : di^Vu ■ VSi = -Hv-Q^ - Q(t) ■ V5i, u = u^ (39) 

(we remind the reader that {VSkl 7^ 0, A; = 1, 2). It will be shown below that the form of the 
vector function Q(t) in f l39|l plays a crucial role. 

First of all, we examine the invariance with respect to the one-parameter Lie groups forming 
the group Gg. Obviously, conditions f l5^ are invariant with respect to the groups Ti,T2 and 
T3 for the arbitrary smooth vector function Q(t) . 

To be invariant under the group To, the conditions must take place 



(99* 

diyVu* ■ VS*i + H,-^ + Q(r) ■ vsi 
at* 



= 0, - w^|(39) = 0, (40) 

(39) 



which lead to the requirement 

Qa{t + So) = Qa{t), a =1,2,3. (41) 
Since fHTj) must hold for arbitrary real values of t and et, we conclude that 

Qa{t) = Qa, a = 1,2,3, (42) 



where qa are arbitrary constants. Thus, the BVP of the form f lT4j) - f[T8|) is invariant with respect 
to the Lie group Tq if and only if conditions ( H2|) take place. 

In a quite similar way, one examines the group T4 and, as a result, arrives at the requirement 

Qaite^'')e'' = Qa{t), a = 1,2,3 
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what implies 



g„(t) = -^, a = 1,2, 3. (43) 



Thus, the BVP of the form ( Tl4|) -( 1T8|) is invariant with respect to the Lie group T4 if and only 
if conditions ( H3|) hold. 

Prior to examine the invariance of the boundary condition (fT6|) with respect to the group 
Tab, we find how the first derivations of the variables u and Si are transformed: 

du* du du du* du du du* du 

^— = Sm Eab + ^— cos Eab, TTT = ^ COS €ab + ^— Sm Eab, TTT = (44) 

OXl OXb OXa OXl OXb OXa OX* OXc 

dSl dSi . dSi dSl dSi dSi . dSl dSi dSl dSi 

= ^ cos Eab, = COS Eab + Eab, ^ = ^> ^ = ^- 

(7X* OXb OXa OXl OXa OX^ OXc Ot* Ot 

(45) 

Substituting (Hlj) and f l45]) into f l4U]) and making relevant calculations, we arrive at the equality 

dSi dSi d Si 

{Qa{t) COSEab + Qb{t) SlYlEab) ^ + *^^^«*^^) siu £„6 + Qb{t) COSEab) 'q^^^''^'^^^ ^ Q(t)-VS'i, 

leading to the algebraic equations 

Qa{t) COSEab + Qb{t) ^VClEab = Qa{t), Qa{t) siu e^b " Qb{t) COSEab = -Qb{t). 

Since Eab is an arbitrary parameter, we immediately conclude that 

Qa{t) = Qb{t) = 0. (46) 

It means that the BVP of the form f|T4|) - f|T8|) is invariant with respect to the Lie group Tab if 
and only if conditions fHBl) hold, while the function Qc{t) is an arbitrary smooth function, i.e. 
the vector function Q(t) contains only one non- vanish component. 

Thus, we proved that the BVP of the form f|T4|) - f|T8|l is invariant under Tq^T^ and Tab iff 
the vector function Q(t) satisfies the restrictions fj42l) . (j43j) and fH6|) . respectively. To finish 
the examination of the BVP with arbitrary functions di{u) and d2{v), we need to investigate 
the case of the arbitrary one-parameter Lie group from the group = (Tq, T4, T21, T31, r32). 
According to the general Lie theory, each one-parameter group of the group ^5 corresponds to 
a linear combination of the generators Pt, Dq and Jab of the form 

L = Ai J21 + A2 J31 + A3J32 + XiPt + Xl + ... + Xl^O, 

where Ai, . . . , A5 are the given real constants. 

It should be stressed that the form of the operator L can be simplified by using the trans- 
formations of variables Xa, a = 1,2,3 from the group E^^^ , namely the linear combination 
A1J21 + A2J31 + A3J32 may be reduced to a single operator of rotation, for example, to the 
operator J12. Indeed, using the transformations 

t^t, X Ai{Pi)A2W2)AM)X 

where the rotation angles a = 1, 2, 3 satisfy the conditions 

Ai cos (32 sin /33— A2 (cos /3i cos — sin (3i sin (32 sin (33)— X3 (sin (3i cos (3^ + cos (3i sin (32 sin (3^) = 0, 
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Ai sin /32 — A2 sin /3i cos /32 + A3 cos /3i cos /?2 = 0. 
one simplifies the operator L to the form 

L = AJ21 + XiPt + X5D0, X^ + XI + XIj^ 0, 

where 

A = Ai cos (32 cos /?3+A2 (cos /?i sin (3^ + sin /3i sin /32 cos /?3)+A3 (sin /?i sin /33 — cos (3i sin /32 cos /33) . 

Now one can easily see that there are only two cases (depending on the parameters A, A4 
and A5 !), leading to new one-parameter invariance groups of the BVP of the form f[T^ - f[TBl) 
with the correctly specified vector function Q(t). This occurs iff: a) AA4 7^ and A5 = 0, and 
b) AA5 7^ 0. 

Let us consider the most troublesome case b). Without loss of generality, one can assume 
A4 = and A5 = 1, therefore, the corresponding one-parameter Lie group is (see Remark 3). 
Substituting transformations from this group into invariance conditions f HU]) . we obtain the 
system 

Qi cos(A£:7) - Q2 sin(A£7) = Qi{t), 
Qi sin(A£:7) + Q2 cosiXej) = Q2{t), 

which can be rewritten as follows: 

Qi = Qi{t) cos(A£7) + Q2{t) sin(A£7), 

Q2 = -Qi(t) sin(A£7) + Q2(t) cos(A£7), 
Q3 = Qsit), 

where Qa = Qa{e'^'^^t)e'^'' , a = 1,2,3. So unknown components of the vector function Q(t) can 
be found from the system of functional equations obtained above: 

qiCOST + q2smT -gi sin r + ^2 cos r ^ _ qs 

Vll^j - , ^2[t) — , ^3[t) — [^1) 

where r = |A logt, qa, a = 1, 2, 3 are arbitrary real constants obeying the condition qiq2X 7^ 0. 

Hence, the BVP of the form f|T^ - f|T8l) is invariant with respect to the Lie group if and 
only if the conditions f H7|) hold. Case b) has been completely investigated. 

Case a) has been examined in a quite similar way. As result, we have proved that the BVP 
of the form fll4p - (IT5]) is invariant with respect to the Lie group Tg iff the conditions 

Qi{t) = Qi cos(At) + q2 sin(At), Q2{t) = -qi sin(At) + q2 cos(At), Qsit) = q^, qiq2X ^ 

take place. 

Thus, we conclude from the analysis carried out above that the trivial group of invariance 
of the BVP of the form f|T^ - f|T8|) is the three-parameter Lie group, generated by the groups 
Ti , T2 and T3 and one is listed in case 1 of Table [2j G'^ is extended either to four-parameter 
group (cases 2-4 of Table [2]) or five-parameter group (cases 5 and 6 of Table [2]) if and only if the 
relevant conditions at the vector function Q(t) take place. These conditions have been derived 
by comparing the components of Q(t) from formulae fH2]) . (1431) and fl46p and those from the 



15 



examination of cases a) and b). One easily notes that there are no other possibihties to extend 
the group G^. It means that the case of the arbitrary functions di{u) and d2{v) is completely 
investigated. 

The next step of the proof is to show that cases 2-9 from Table [T] do not lead to any 
new symmetries of the BVP of the form (IT^ - (fT5]) . Consider case 2 when di{u) = ki and 
d2{v) is an arbitrary function. According to Theorem 1, the NHE system (IT^ - (IT^ admits the 
infinite-dimensional Lie algebra formed by the operators from the basic algebra AE{1,3) and 
the operators udu, a(t, x)5„. The corresponding groups to these operators are generated by the 
transformations 

T[/ : t* = t, X* = Xa,a = 1, 2, 3, u* = e^^u, v* = v, = Si, 5*2 = 5*2, 

Too : t* = t, X* = Xa,a = 1,2,3, u* = u + a{t,x.)eoo, v* = v, SI = Si, S2 = S2, 

Let us prove that any BVP in question is not invariant with respect to the group Tu. In fact, 
the second invariance condition in (HOj) takes the form 

and can be satisfied for an arbitrary value of the parameter eu only in the case = 0. On the 
other hand, if = then Um 7^ and, thereby, the second boundary condition in f|T7|l is not 
invariant with respect to Tu. In a similar manner, it can be shown that any BVP in question 
is not invariant with respect to Too. 

Now, we must examine each one-parameter group, corresponding to a linear combination 
of the operators Pt, Dq, Jab, b < a = 1,2,3, udu and a{t, x)(9u (see case 2 of Table [1]): 

Y = Ai J21 + A2 J31 + A3 J32 + XiPt + X5D0 + \&udu + \7a{t, x)9„, A^ + A^ ^ 0, 

where Ai, . . . , A7 are arbitrary real constants. To avoid cumbersome formulae, we present in an 
explicit form only the point transformations for the variable u: 

u* = e^^'^u + A7 /" a*{T)e^'^^-'^^dT, (48) 
Jo 

where a*(r) = a (r (r), x*(r)). Obviously, if Ae = 0, A7 7^ or As 7^ 0, A7 = 0, then the BVP 
under study is not invariant with respect to the relevant one-parameter Lie group because one 
uses the result obtained above. Hence, we should examine only the case A6A7 7^ 0. Taking into 
account formula (jlSl) . the second invariance condition from fHU]) takes the form 

A7 / a*(r)e^«("-^^)dr = {l - e^«"^) . 
Jo 

Thus, we obtain 

Using this formula, the invariance criterion for the second boundary condition from f[T7|) leads 
to the condition 

(1 - e^«^^) =u^{l~ e^^'^) . 

Since this equality must hold for arbitrary values of the group parameter Ey, we immediately 
obtain that Um = u^. It is nothing other than the contradiction because Um 7^ and Ag 7^ 0. 
Thus, we conclude that the operator Y is not a Lie symmetry operator for any BVP from class 
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f|T^ -f lT8|) . Case 2 from Table |2] is completely examined. Obviously, case 3 from Tabled] can be 
studied in a similar manner to case 2 (the boundary conditions f|T7|l and f|T8l) should be used) 
and the same result is obtained. 

Consider cases 4 — 8 from Tabled] when di{u) and d2{v) are specified non-constant functions. 
It turns out that cases 4 — 7 from Table d] can be examined in a similar manner as we did it in 
case 2 using the group Tu. Finally, no new Lie groups of invariance are obtained. 

The most non-trivial is the case of conformal power n = m = (see case 8 in Table dj) 
because of the conformal operators Kb = \:K\'^dx^ — 2xbXadx^ + 5xbudu + bxbvd^, 6 < a = 1, 2, 3, 
which generate the one-parameter Lie groups 

rp _ * Xa ^ Xh |X| Sfs-^ ^ X, 

->- Kh ■ X„ — - ; TTT-;^ , X 



Xb - 


X 




1 - IXbEK^ + 


X 





1 - IXbSKb + ' 1 - IXbEKb + ' " 1 - IXbEK^ + jx^e^^^ ' 

5 5 

f = t, u* = u{\- 2xbeKi + Ixp^xJ ' > V* = v{l- 2xbeKi + Ixl^e^J ' , 

where a,b, c = 1,2,3; a ^ h,a ^ c,h ^ c. One notes that the boundary condition ( 1T8|) is not 
invariant under T^^^. In fact, the invariance conditions 

lim |x*| = -Fcx), 'y-Woo||m = 0. (49) 

|X|^+oo 

are not satisfied because 

lim |x*| = lim ^ — ^ = +oo. 

ixK+oo ixK+oo _ 2^^^^^ + |x|%2^J ^ 

Thus,we conclude that any BVP of the form f|T4l) - f|T8l) is not invariant with respect to the Lie 
group Tk^. 

Finally, we have examined all possible one-parameter groups, corresponding to linear com- 
binations of the operators Kb, b = 1, 2, 3 and D (see case 8 in Table d]) 

Z = XiKi + X2K2 + X3K3 + X^D, Xl + Xl + Xjj^ 0, 

and shown that the boundary condition (1181) is not invariant under those groups. 

Nevertheless, the case of linear governing system fll4l) - (ll5l) (case 9 in Table has a cum- 
bersome algebra of invariance; it was examined in a quite similar way to case 2 and no new 
invariance groups found. 

The proof is now complete. ■ 

Remark 4 In Theorem 2, the special case of system [14\)- [T5\) with di{u) = d2{y) (see case 10 
in TableUD was not examined because the relevant BVP with the equal dijfusivities is unrealistic 
from the physical point of view. In fact, the vast majority of substances have different physical 
characteristics of solid, liquid, and gas phases. 



4 Symmetry reduction and invariant solutions of BVPs 
from class (1141)— (flSl) with the constant energy flux 



4.1 Optimal system of subalgebras of the invariance algebra 

Let us consider a nonlinear model of heat transfer processes in metals under the action of 
intense constant energy flaxes directed perpendicular to the metal surface. This model coincides 
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with the BVP ©-([n]) (see also (HD-dlH])), where Q(t) = q = (0,0, g), q = const. According 
to Theorem 2 such a BVP admits the five-parameter Lie group G5 of invariance (see case 5, in 
Table E]) formed by the one-parameter groups Tq, Ti, T2, T3 and T12. This group corresponds to 
the five-dimensional Lie algebra with the basic operators 

Pt = Pa = ^— , Ji2 = X2P1 - X1P2, a = 1, 2, 3. 

at OXa 



Our primary aim is to show how using these operators one can reduce the BVP (fT^ - (fT8|) . 
where Q(t) = q, to BVPs of lower dimensions. We use for this purpose the optimal systems 
of s-dimensional subalgebras {s < 5) of A^. All such subalgebras are non-conjugate up to 
the group of inner automorphisms of the group G5. To construct a full list of the optimal 
systems, we represent the algebra A5 as follows: A^ = {Pi, P2, J12) ® {P3) © (Pt)- Now by using 
the well-known Lie-Goursat classification method for the subalgebras of algebraic sums of Lie 
algebras [3^ (see the monograph |1S] for details) and the results of the subalgebras classification 
of low- dimensional real Lie algebras [IH] , it is easy to obtain the complete list of subalgebras of 
the algebra A^. This list can be divided into subalgebras of different dimensionality. 

One-dimensional subalgebras: 

(P3 cos + Pt sin 0) , {Pi + a (P3 cos (p + Pt sin (p)) , {Ju + /3 (P3 cos (p + Pt sin 0)) ; 

Two-dimensional subalgebras: 

(P3, Pt) , (Pi + a (P3 cos (p + Pt sin cp) , P2) , (Pi + a (P3 cos ^ + Pt sin 0) , P3 sin (p - Pt cos (p) , 

( J12 + /3 (P3 cos (p + Pt sin (p) , P3 sin (p - Pt cos (p) ; 
Three-dimensional subalgebras: 

(Pi, P3, Pt) , ( J12, P3, Pt) , (Pi + « (P3 cos <P + Pt sin (P) , P2, P3 sin - P^ cos 0) , 

( J12 + (3 (P3 cos (P + Pt sin 0) , Pi, P2) ; 
Four-dimensional subalgebras : 

(Pi, P2, P3, Pt) , (J12 + (P3 cos (P + Pt sin 0) , Pi, P2, P3 sin - P^ cos 0) ; 

Five- dimensional subalgebra: 

(•^12, Pi, P2, P3) Pt) ■ 

where a > and /3 are arbitrary real constants, < < tt. 

We remaind the reader that the additional conditions on the functions Sk{t,x.), Yk{t,:K) 
and Q(t) arising in the BVP class f lT^ -f lT8|) have been imposed, namely 

dS 

^ ^ 0, \VSk\ + 0, Q(t) ■ m ^ 0, Vfe ■ rifc ^ 0, = 1, 2. (50) 

On the other hand, the Lie-Goursat classification method is a purely algebraic procedure; 
therefore, some subalgebras presented above lead to invariant solutions, which do not satisfy 
the restrictions fl50|) . For example, the algebra (Ji2,P3,Pt) generates the ansatz 

u = u{r), V = v{r), Sk = Sk{r), k = 1,2, 
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where r = \/x'l + x^. Obviously, one sees that ^ = and q ■ ni = 0. Thus, the contradiction 
is obtained and we conclude that the algebra {Ju, P3, Pt) leads to solutions, which do not have 
any physical meaning. 

The complete list of the subalgebras, leading to invariant solutions of the BVPs in question 
satisfying the restrictions f l5U]) . reads as follows. 
One-dimensional subalgebras: 

(P3 cos (t) + Pt sin 0) [j) 7^ 0, 1) , (Pi + a (P3 cos + P* sin 0)) , ( J12 + /? (P3 cos + P* sin 0)) ; 

Two-dimensional subalgebras: 

(Pi + a (P3 cos + Pt sin 0) , P2) , (Pi + a (P3 cos + P* sin 0) , P3 sin - P^ cos 0) [j) 7^ 0, ^ , 

(Ji2 + /3(P3Cos0 + Ptsin0),P3sin0-Ptcos0) (^0 7^0,^) ; 
Three-dimensional subalgebras: 

(Pi + a (P3 cos + Pt sin 0) , P2, P3 sin - P* cos 0) ^0 ^0,0, (51) 

{J12 + P (P3 cos + Pt sin 0) , Pi, P2) (0 ^ 0, 1 if /? ^ 0) ; (52) 
Four- dimensional subalgebras : 

( J12, Pi, P2, P3 sin - Pt cos 0) (0 7^ 0, ^) , (53) 

where a > and (5 are arbitrary real constants, < < vr. 



4.2 Symmetry reduction and invariant solutions 

Now, one may use each algebra from this list for reducing the BVP of the form f[T^ - f[T51) 
with Q(t) = q to the BVP of lower dimensionality and to construct the exact solutions of the 
problem obtained. We also note that three- and four- dimensional subalgebras (l52ll and (l53l) 
generate the same ansatz to find the unknown functions m, Si and 5*2: 

u = u{z), V = v{z), Sk = Sk{z), k = 1,2, z = X3 — fit. (54) 

(Hereafter, we use the designation fi = — tan0). The application of three-dimensional subal- 
gebra (l5Ti) also yields ansatz (15^ but with the invariant variable z = a*xi + x^ — fit, a* e R. 
Ansatz ([51]) reduces the BVP of the form ([Il])-([l8]) to the BVP for ODEs, which was studied 
in details in our earlier papers [22l[23l[35] . This ansatz leads to the plane wave solutions of the 
BVP in question, while melting and evaporation surfaces are two parallel planes moving with 
unknown velocity fi along the axes Ox^. 

New non-trivial reductions occur if one applies one- and two-dimensional subalgebras. Let us 
consider the algebra ( J12 + /? (P3 cos + Pt sin 0) , P3 sin — Pt cos 0). Solving the correspond- 
ing system of characteristic equations, one obtains the ansatz 

u = u{r, z), V = v{r, z), Sk = Sk{r, z), k = 1,2, (55) 

where z = x^ — fit — /3 arctan f^, = ^/ xf + X2 are the invariant variables. 
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Substituting ansatz fl55|l into the BVP of the form f|T^ - f|T8|l with Q(t) = q and making the 
relevant calculations, we arrive at the two-dimensional BVP 



dSi 

Si{r, z) = : diyVu ■ VSi = {nHy - q) u = u^, (58) 

882 

S2{r, Z) =0 : d2mVv ■ V S2 = dimVu ■ V'S'2 + flHm-^, U = Um, V = Vm, (59) 

r'^ + z'^ = +00 : V = Voo, (60) 



where is an unknown parameter and V'= l^^iyTT + l^ 

Nevertheless, the BVP of the form fl56p - (H7]) is much simple object than the initial BVP, 
but it is still the nonlinear problem with the basic two-dimensional PDEs. Our purpose is to 
reduce one to a BVP with basic ODEs. Of course, one may apply different technics to realize 
such reduction; however, we confine ourselves to the case when the invariant variables r and z 
admit clear physical meaning. It happens for (3 = because then the first variable z makes the 
transition to a moving coordinate system (in the direction of the variable X3) with the origin 
at the evaporation surface, while the variable r presents the radial symmetry of the process 
with respect to the variables xi and X2- Obviously, such a situation takes place if the surface 
bounded by a circle of the radius R is exposed by the flux Q(t) = q. 

Thus, setting /? = 0, we may consider the ansatz 



u 



u(u), V = v{uj), Sk = Sk{co), u> = z + V z^ + r^, k = 1,2, (61) 



used earlier in [31] for a similar purposes. Note that it is a non-Lie ansatz because the maximal 
algebra of invariance of system f l56p - fl57l) (with arbitrary non-constant functions di{u) and 
d2{v)) is trivial and generated by the operator 

Remark 5 Recently we found that ansatz with oj determined from the cubic equation 

l--(x3-/it) = ^+ , . ,/teR (62) 

was used in paper 15^ to construct the exact solution of the BVP of the form [l4\)- fT^) with 
Q{t) = q and the constant difjusivities di{u) and d2{v). Setting k = and — fit = z in [6^] 
one arrives at ansatz (EIP- However, we have checked that ansatz with u defined from 
[6^) ( with any non-zero k, !) is not applicable to reduce the BVP of the form [l4\)- fT^) with 
non-constant diffusivity. 



Substituting ansatz fl6T|l into the BVP of the form fl56|) - fl60|l and taking into account the 
relations 

^ 2uj du dSk ^ 2u dv dSk dSk uj dSk , , „ 

Vm ■ V'^fc = Vv ■ VSk = r^, ^ = — ^, k = l,2, 

\/ Z^ _|_ ^^2 _|_ J.2 du OZ ^J + du 
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we obtain the BVP for ODEs: 

d f . . ,du\ u du ^ ^ 

— \ujdi[u)— ] + II-— = 0<wi<w<a;2, (63) 
auj \ du J I du 

d f ^ , ^ dv\ u dv , ^ 

du 

u = Ui : 2diy— = uHv — q, u = m„, (65) 
duj 

dv du 

uj = UJ2: 2d2m-r- = 2rfim— + nHm, U = Um, V = Vm, (66) 
du du 

u = +00 : V = Voo, (67) 

where Uk, k = 1,2 and fi are parameters to be determined. 

Now, we can define the form of the free surfaces Sk{t, x) = 0, k = 1,2 because in accordance 
with ansatz fl6T]) 



Sk{u) = z + V z"^ + r"^ = Uk, k = 1,2. 
Obviously, the last equations can be rewritten in the form 

^ , ' = 1 , k=l,2. 68 

Thus, the equations obtained define some paraboloids of revolution in the space of variables 
Xi,X2,z (see Fig.l). From the physical point of view, unknown parameters should satisfy the 
inequalities > > 0. Moreover, the parameter ui can be defined as follows. If one sets 
2; = in fl68l) then ui = \/xl + x\. On the other hand, a part of the surface bounded by a 
circle of the radius R is only exposed by the fiux Q(t) = q, i.e. we can set ui = R without 
losing the generality. 

Let us now turn to the construction of the exact solutions of problem (jSSD^dSZD- In fact, 
the general solution of the nonlinear equations (163 p and fl64p is unknown. However, in some 
cases, one is known (see, e.g., [51]). Here, we consider two cases in details. 

Example 1. The BVP of the form ( fg3j) -([g7|) with di{u) = ai and d2{v) = a2, where 
ai , as e M+ . 

In this case, the general solutions of equations flB51) and fl64p are given in an explicit form 
by the formulae 

u = Ci^iiu) + C2, v = Cs^2iuj) + C4, (69) 

-[-00 U) 

where $fc(a;) = J^°° u^^e du, k = 1,2, Ci,...,C4 are arbitrary to-be-determined con- 
stants. Substituting solution (169|) into the boundary conditions (l65l) - (l67|) and taking into 
account the formulae 

— — = u e , k = 1, 2, 
du 

we arrive at the exact solution 

Uy-Um ^ / N , Um^i{R) - U^^i{u2) Vm - V^o , . , 

" = ^ /m ^77 — H ^ -7—, — ^ — , V = —— — ^$2(w) + ^;oo, (70) 

$l(i?) - $1(^2) ^i{R)~^i{u2) $2(^2) 

where the parameters U2 and must be found from the transcendent equation system 

^^1" ^ /m — -i?~^e"^^ = fiH^ - q, 

^i{R)-^i{u2) 

^2{i02) ^i{R)-^i{u2) 
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Finally, using formulae fl70|l , f l55|l and fl6T|) . we obtain the exact solution of the BVP of the 
form f ll4l) - f|T8l) with di{u) = ai, d2{v) = 02 and Q(t) = q in the explicit form: 



V 



$2(^2) 



$2 x\+xl + (Xs - + - llt^ + ■Uoo, 

x\ + xl 2{x-i — lit) 
Sk = ^ , ^ + ^ -1 = 0, /c = 1, 2; = i?. 

Remark 6 Whereas the basic equations of the BVP of the form [l4\)- fT^) with di{u) = ai, 
d2{v) = 02 are linear, the relevant equations of the initial BVP of the form ^-IjJ^ may be 

nonlinear, but satisfvinq the conditions ai = tttt^tj-t^, 0.2 = ^^1*^-/7!! [see the functions (hi and 
' jy ^ Ci(<^/(«))' ^ C2(02 (")) ^ 

02 m (d;. 

Example 2. The BVP of the form l6^) -^6l\) with di{u) = u^^ and d2{v) = 1, i.e., (E^j is 
the fast diffusion equation, while [64\ ) is the linear diffusion equation. 

In this case, the general solutions of equations ( 163!) and (l64l) are given by the formulae [51] 



Z""^" du 

I ^7(YT^r^^^ = in- + C2, . = C3$M + Q, (71) 

where $(aj) = J^°° u^^e^^'^du, W{x) is the Lambert function, A = — |z/ + Ci, a is an arbitrary 
constant, Ci, . . . , C4 are to-be-determined constants. 

Substituting solution f lTT]) into the boundary conditions flB^ - flB7|) and taking into account 
the formulae 

(i$ 1 fj.,, ^ /udu\ u du 

- u~^e~^'^, In --— + --— = A, 



du \u du J u du 

we obtain the exact solution of the BVP in question 



In— , w = — - — —^{u) + v^, (72) 



(1 + e-^('=^)+^) i?' ^U2) 
where the parameters 072 and /i must be found from the system of transcendent equations 



^{U2) 



Here, we used the designations 



A = + In I [^H^ ~^'^) ^ " ~ 2" ^ 2 
A{u2) = --uj2Um + In I {iiH^ - g) 17 I + {^^H^ - g) 77 + o^^''- 
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Finally, using formulae fl72|l , f l55|l and fl6T]) . we obtain the exact solution of the origin BVP 
f lT^ -f lT8|) with di{u) = u"^, d2{v) = 1 and Q(t) = q in the implicit form 



It should be noted that several BVPs of the form (l63|) - (l67|) can also be exactly solved 
for some other forms of diffusivity coefficients. Moreover, one may apply the standard pro- 
gram package (e.g.. Maple, Mathematica) to numerically solve this BVP with arbitrary given 
diffusivities. 

5 Conclusions 

In this paper, multi-dimensional nonlinear BVPs with the basic evolution equations by 
means of the classical Lie symmetry method are studied. We consider BVPs of the most 
general form , which include basic equations of the arbitrary order, boundary conditions 

on known and unknown moving surfaces, boundary conditions on regular and non-regular 
manifolds. A new definition of invariance in the Lie sense for such BVPs is formulated. The 
definition generalizes those proposed earlier for simpler BVPs [2l[33l|35l|36], and it can be 
extended to BVPs for hyperbolic and elliptic equations. Note that the comparison of this 
definition with those proposed earlier is presented in the recent paper [53], where Definition 1 
in the case of (1+1) -dimensional BVP was formulated. In this paper, we also propose the 
algorithm of the group classification for classes of BVPs, i.e., extending the well-known problem 
for PDEs to BVPs. Of course, the group classification problem for simple classes of BVPs can 
be solved in a straightforward way (see, e.g., |24j), however, one needs to determining some 
algorithm in the general case. 

The main part of the paper is devoted to solving the group classification problem for the 
class of (l+3)-dimensional BVPs fl8|)- f|T2|) . modeling processes of melting and evaporation under 
a powerful fiux of energy. First of all, we simplified the BVPs in question to the form f lT^ - f|T8l) 
using the Goodman substitution. Since the basic equations of the problem obtained are the 
standard NHEs, we used the known system of determining equations [111110] to derive their Lie 
symmetry description. Having done this and using the group of equivalence transformations, 
we proved Theorem 2 presenting all possible Lie groups of invariance of the BVPs of the form 
( IT4|) -( [T8|) . It was shown that the Lie invariance does not depend on the diffusivities di{u), 
d2{v) but only in the form of the flux Q(t). There are only flve correctly specifled forms of 
Cl{t) (see Table 2) leading to the extensions of the three-dimensional Lie group of invariance 
(the trivial Lie group), which is admitted by the arbitrary BVP of the form (I14l) - fll8l) . One 
may note that cases 3 and 4 from Table 2 have no analogs in the (l+l)-dimensional space of 
independent variables, while cases 5 and 6 are generalizations of the relevant (l+l)-dimensional 
BVPs (see [24J for comparison). 

We study in detail the BVP of the form f[T^ - f[TSl) with arbitrary diffusivities and the special 
form of the flux Q(t) = q, which naturally arises as a mathematical model of the melting and 




^yxl + xl + {Xs - IJ,ty + Xs - fit 

R 



OO ! 
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evaporation process. Since the MGI of this problem is five-dimensional, the sets of optimal s- 
dimensional subalgebras were constructed using the known algorithm [171H8] . The brief analysis 
of these subalgebras is presented. One of them, the two-dimensional algebra, is applied for the 
reduction of the problem in question to the nonlinear B VP for ODEs. Finally, the BVP obtained 
was exactly solved in two cases of correctly-specified diffusivities; hence, the exact solutions of 
the BVP of the form (fT^ - (|T8|) with these diffusivities were found. It should be noted that 
the BVP of the form f|T^ - (|T8l) with constant diffusivities was studied earlier in |50], where 
the same result was obtained using an ad hoc ansatz, which does not connected with any Lie 
symmetry. To the best of our knowledge, the exact solution of the BVP of the form (fT4l) - (fT8ll 
with the fast diffusion constructed above is new. 

The work is in progress to apply the results obtained in this paper to the reduction and 
construction of exact solutions for other multi-dimensional BVPs with the remarkable Lie in- 
variance. 
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